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Abstract. A mass correlation of central black holes and 
their spheroids ~ 0.002 (within a factor of three) is sug- 
gested by Hubble Space Telescope (HST) and various 
ground-based CCD photometries of early type galaxies. 
The near-IR images of quasar hosts and the emission 
line measurements of Broad Line Region for bright QSOs 
present a similar correlation, which supports the specula- 
tion of an evolutionary linkage between the early active 
QSO phase and the central black holes in normal galax- 
ies. On the other hand, recent reverberation mapping of 
a sample of Seyferts shows a broad distribution of black 
hole to bulge mass ratio with a mean of ^ lO^"^'^, about 
one magnitude lower than the value in early type galax- 
ies and bright QSOs. Adopting a simple unified formation 
scheme for QSOs and Seyferts, we will discuss in this let- 
ter the dependence of the black hole to bulge mass ratio 
in Active Galactic Nuclei (AGNs) on the environmental 
parameters of the host galaxies. We show a broad distri- 
bution of the mass correlation could be due to different 
velocity dispersion of the accreting gas from different for- 
mation mechanism, and the mass ratio in normal galaxies 
and bright QSOs is probably a limit case of black hole evo- 
lution by merger enhanced accretion close to Eddington 
limit. 
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1. Introduction 

The existence of supermassive black holes, at least in ellip- 
tical and bulge-dominated galaxies is suggested by various 
observations, such as the optical spectroscopy, VLBI wa- 
ter maser measurements and X-ray observations of AGNs 
(Nakai et al. 1993, Kormendy & Richstone 1995, Faber et 
al. 1997, Miyoshi et al. 1995, GreenhiU et al. 1996, Pounds 
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et al. 1990). An intensive discussion about the possible re- 
lationship between the central properties and their host 
galaxies was spurred on especially since the high resolu- 
tion HST photometry and several ground-based CCD pho- 
tometries of early type galaxies . Although there are sig- 
nificant scatters, a linear mass correlation of central black 
holes and their host spheroids has still a high probability. 

The theoretical interpretation for such linear scale is 
discussed by several authors (Mcrritt 1998, Silk & Rees 
1998, Wang & Biermann 1998). Considering various obser- 
vations, Wang & Biermann (1998) proposed a possible for- 
mation scenario for elliptical and bulge-dominated galax- 
ies, as well as the consequential active evolution phase 
where early type galaxies are the products of major merg- 
ers between two comparable disk galaxies; the violent col- 
lision between two galaxies could destroy the original stel- 
lar disks and form the spheroidal component of merging 
galaxies after the relaxation; help to release the angular 
moment of the cold gas outside and drive them inwards; 
a central starburst and QSO accretion in the central con- 
densed gas disk will compete for the gas supply, feedback 
and drain the gas in the disk in a short time; possibly 
grow a supermassive black hole during the spheroid for- 
mation, blow up the gas left probably by the nuclear wind 
when the central engine gets to be powerful enough; thus 
restrict the central black hole mass and the mass of the 
stellar component to a ratio of ^ 10^^. 

The numerical simulation by Wang & Biermann (1998) 
shows the star formation approximately scales with the 
nuclear accretion during galaxy interactions, which can 
regulate the black hole to bulge mass ratio to the observed 
level ~ 0.002 within a factor of three. In this scenario, a 
quasar black hole is possibly formed at a cosmic time scale 
of ~ lO^yr, corresponding to z = 3 5 in a flat Einstein- 
de Sitter world model with il = 1, Hq = 50knis^^ Mpc~^. 
Afterwards, when the quasar active phase has ceased and 
the spheroidal stellar system gets relaxed in a comparable 
relaxation time scale, an elliptical or spiral bulge harbor- 
ing a massive black hole in the center with the mass in 



2 



Wang et al. : Black hole to bulge mass correlation 



a factor ~ 10 ^ of the host spheroid may appear in the 
universe. 

We should mention the black hole evolution in this 
model is assumed to reach the Eddington accretion rate 
whenever there is sufficient gas supplied to the center. The 
question of whether this mass correlation is universal for 
all AGNs, i.e. whether Seyfert galaxies follow a similar 
black hole to bulge mass correlation as in early type galax- 
ies and luminous QSOs, is not only important for the un- 
derstanding of the correlation between the host galaxies 
and the evolution of active nuclei, but for the relation of 
the formation and evolution scheme between Seyferts and 
QSOs. 

Recent reverberation mapping of a sample of Seyferts 
with reliable central masses and the bulge magnitudes by 
Wandel (1999), suggests that there is a broad distribu- 
tion of black hole to bulge mass ratio with a mean of 
~ lO^"''^, about one magnitude lower than the value in 
early type galaxies and bright QSOs. Although this dis- 
crepancy may in part be explained by systematic errors 
and selection or orientation effects (Wandel 2000), a sig- 
nificant difference probably remains, in particular between 
Seyfert 1 galaxies and QSOs, at least a large dispersion for 
the medium-bulge system, with masses similar to those of 
Seyferts (Magorrian et al. 1998, Richstone et al. 1998, Fer- 
rarese & Merritt 2000). Our purpose is to firstly present 
a model which could explain the black hole to bulge mass 
correlation in AGNs and the dependence on the environ- 
mental parameters of the host galaxies, such as the gas 
or stellar velocity dispersion, as well as the relation of the 
central starburst and accretion process during galactic in- 
teraction; Secondly, discuss the dispersion of black hole 
to bulge mass ratio in QSOs and Seyferts within a simple 
unified formation scheme, where the bulge formation and 
nuclear activity are triggered by galaxy mergers or tidal 
interactions. We found the variation of the velocity disper- 
sion of accreting gas could cause a range of distribution 
for the mass ratio, leading to a correlation of the nuclear 
black hole mass and the gas velocity dispersion roughly 
as Mbh oc Ug'^ from our simulation, within the slope sug- 
gested by recent work of Ferrarese & Merritt (2000) and 
Gebhardt et al. (2000). 



1995, Taniguchi 1999). In this sense, QSOs and Seyferts 
are the consequences of either galaxy mergers or tidal in- 
teractions, meanwhile a massive black hole may grow dur- 
ing the process of mergers or tidal interactions since the 
enhanced cloud-cloud collision or star formation could in- 
crease the mass inflow to the center. 

We adopt such a unified formation scheme for Seyferts 
where we think the tidal perturbation of some satellites 
not only kinematically heat and thicken the host disks, 
form the central bulge, but may increase the effect of self- 
gravity in the molecular disk, enhance the cloud-cloud col- 
lision and star formation, thus result in a considerable 
mass inflow to the central region on a relatively short 
time (Lin et al. 1988, Quinn et al. 1993, Walker et al. 

1996, Velaquez & White 1999). Although we consider in 
our model the tidal interactions as a formation scheme for 
most Seyferts, we do not reject a possibility that some of 
AGNs with low black hole to bulge mass ratio could still 
be major mergers, but at the earlier evolutionary phase 
in which the black holes have not had time to reach the 
asymptotic value. 

We adopt E{R) = as the initial mass distribution 

for protogalaxies in our model, with the typical mass of 
10^^ Mq and the disk scale of ~ 14kpc, similar as the size 
of our Galaxy. The evolution of the surface density S(i?, t) 
of a differentially rotating disk with angular velocity and 
viscosity u is governed by (Liist 1952, Pringle 1981) 



5E 
'dt 




dR 



d_ 
dR V 



uER' 



dR.) 



(1) 



where u — (3iv^r for a Keplerian selfgravitating disk 
(Duschl et al. 2000); the accretion time scale tacc jv = 
r/Piv^; and the star formation time scale = atacc = 
P2r/v^ (Pringle 1981). S/t* is the star formation rate in 
the disk with mass return rate Rg ~ 0.3 (Tinsley 1974). 
In this case, the mass influx at the inner boundary i?in is: 
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2. Black hole to bulge mass correlation in AGNs 
2.1. Illustrative model 

AGNs are usually thought to be powered by gas accretion 
to the central massive black holes. There; is now firm evi- 
dence that QSOs could be formed during the violent galac- 
tic interactions together with central starbursts (Sanders 
et al. 1988a, 1998b, Taniguchi & Shioya 1998, Heckman 
et al. 1986, BahcaU et al. 1997). Although only ~ 10% 
of Seyfert galaxies have companions, there are still many 
aspects of observations which support a simple unified for- 
mation scheme for QSOs and Seyferts (Rafanelli et al. 



We adopt Vrj, oc Vq for a flat rotation law for the proto- 
galaxies in our calculation with Vq = lOOkm/s as a nor- 
malization velocity. The justification of such simplification 
for a self-gravitating disk at parsec to kiloparsec region is 
discussed in detail by Wang & Biermann (1998). 

The spheroidal Bondi accretion rate is given by m = 



'InXpRi 



where Rare is the Bondi accretion radius de- 



fined as i?acc = (-^bh is the mass of the central black 

hole; Vc the effective relative velocity between seed black 
hole and the ambient gas) (Bondi 1952). In our model, 
we start from a tiny black hole ~ 5Mq, with the inner 
boundary of the accretion disk ~ O.lpc which is usu- 
ally thought to be the inner edge of a torus. Changing the 
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scale of inner boundary or the seed black hole mass only 
influences the early evolution, but no significant effect on 
the final result. At early time, Rm » Race, we assume the 
mass shears inwards at the inner radius to form a uniform 
Bondi flow with the mass distribution = ^""^^3°^^^° 

{At is the time step of the calculation); afterwards, black 
hole would grow and reach a stage of Race ^ -^in; we as- 
sume in this case black hole accretes the mass shearing 
inwards within its influence at i?acc via a uniform Bondi 
flow with the mass distribution = ^4^3^°^'^° ■ 

The Bondi parameter is assumed to be A = 1/2 in 
our calculation as a possible reduction factor due to the 
angular momentum. Thus the accretion rate rh and the 
star formation rate 'J in a time step At is given by: 

3 A /?! 2 ^ So Ro M^jRl v^, if i?acc < Rm 
m= ^ (3) 
3 A /3i 1/0 2 ^ So i?o v^/G Mbh if i?acc > Rin 
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We could roughly estimate the black hole to bulge mass 
ratio by f = GAf.^wrf/fiij ^ ^^ere we take the value 
of m in case of i?acc > -Rin- The estimated ratio shows a 
strong dependence on the velocity dispersion Uq, the value 
(3i * /32 = t*/t acc and the black hole mass Mbh- 

From Eq. (^), we know the accretion rate to oc 1/Mbh 
when -Race > ^^in, which indicates that the nuclear activity 
would dim slowly as the central black hole becomes mas- 
sive enough. In this case, a critical black hole mass would 
be reached when the nuclear accretion rate decreases to a 
level below Eddington rate medd — 2.16 lO"'' [M^h/MQ) 
in units of to ^ lO^yr. Assuming to ~ TOodd, we get the 
critical black hole mass ^ 2.6 IO'^Mq, with Vc ^ 50km/s 
and (3i ^ 0.05 corresponding to the viscous diffuse time 
scale of ~ lO^yr (we will discuss these particular parame- 
ters later). Afterwards, the accretion rate would decrease 
quickly to be much less than Eddington rate, and the 
black hole will not grow significantly. We say this is an 
upper limit, simply because we use the initial surface gas 
density S(_R) = throughout the estimation. Actu- 

ally, the gas would be depleted during the evolution by 
star formation and accretion etc., it would thus take a 
time scale much longer than Hubble time in order to ap- 
proch such a mass or maybe never could reach such a 
level if there is no external trigger for the nuclear accre- 
tion, such as galaxy mergers or tidal effects. In any sense, 
the black hole growth would be locked up in a level of 
~ 10'' Mq in such an evolutionary scheme. Meanwhile, we 
assume in our model /3i/?2 ^ 1, since the result of nu- 
merical simulation by Wang & Biermann (1998) shows a 
strong correlation between the starburst and accretion in 



merging galaxies, which could regulate the black hole to 
bulge mass ratio to a level ^ 0.002 within a factor of three 
in QSOs and early type galaxies. In this case, if we con- 
sider Seyferts and QSOs are formed in a similar scheme, 
we could adopt this result as a reasonable assumption for 
our present model. We now use as input the critical black 
hole mass 2.6 10^ M©, the assumption /3i /?2 ~ 1 and 
Vc ^ 50km/s to Eq. (||), a mass ratio about 0.001 would 
be approched. Choosing Vc ~ 50km/s is based on the ob- 
servations by Zylka et al. (1990) of molecular line emission 
in the central region of our Galaxy with the assumption 
that the velocity field of molecular clouds would be the 
low limit of the sound speed of hot gas (Linden et al. 1993, 
Garcia- Munoz et al. 1977, Reynolds 1989, 1990). We know 
from Eq. (^ and the numerical calculation, the variation 
of the velocity dispersion of accreting gas, Wc, would cause 
a broad distribution of black hole to bulge mass ratio, 
which is actually shown by the observations. The ratio 
could reach ^ 0.06 when ^ 180km/s, close to the virial 
velocity of the system. In this case, the black hole evolu- 
tion would follow Eddington accretion for a longer time, 
resulting in a QSO phase finally. We should mention here 
the range (0.001 ~ 0.06) is a crude upper limit, since we 
could not include the surface density evolution by star 
formation and accretion in this estimation. The practical 
values should be given by numerical simulation shown in 
Fig. I 

2.2. Numerical results and discussion 

To numerically solve the partial differential Eq. ([^) , we in- 
troduce the dimcnsionless variables, r' — r/r^^ t' = t/to, 
S' = E/So with the scale ro = Ipc, to - lO^'yr. The 
general boundary conditions /iTOr_»i?,,„jS(r, t)] = 0, and 
zero torque at the origin, lirrir^B,^^ [G{r^ t)] = 0. The vis- 
cous torque G is given by Pringle (1981) as G{r,t) = 
2TTiyT.R^^, with the disk's outer radius Rout ^ 14kpc 
and the inner radius R^ ^ 0.1 pc. 

This choice of boundary conditions guarantees zero vis- 
cous coupling between the disk and the central object, al- 
lows all mass reaching the inner boundary to flow freely 
inward. 

The numerical results are shown in Fig. 0, where we 
model the effects of a tidal perturbation in a few kpc re- 
gion by increasing the viscous friction. In fact, the numer- 
ical simulation by Lin et al. (1988) demonstrated that a 
tidal disturbance could propagate to the nucleus from a 
distance of a few kpc on a time scale shorter than the 
tidal interaction by "swing amplication" , which induces a 
rapid increase in the effective viscosity and results in an in- 
crease in the rate of mass transfer and energy dissipation. 
We adopted in the calculation a reasonable set of param- 
eters which could match the observations quite well. The 
viscous diffuse time scale is ~ lO^yr (/3i ^ 0.05) in case of 
a tidal perburbation according to the numerical result of 
Lin et al. (1988). We consider the tidal perturbation would 
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Fig. 1. The black hole to bulge mass ratio versus velocity 
dispersion of accreting gas to the central region from our 
simulation. The square corresponds to a mass ratio at the 
level of the mean value of Seyfert 1 samples from reverber- 
ation mapping by Wandel (1999) and the velocity disper- 
sion of ^ 50 km/s; The star corresponds to a mass ratio 
~ 0.002, about the observed value in QSOs and normal 
galaxies, with Vc ~ 180 km/s close to the virial velocity 
in the bulge system in our calculation. We could obtain 
a rough proportionality of the black hole to bulge mass 
ratio versus the veloctiy dispersion of the accreting gas as 
Mbh/Mbuigc oc v^-^. 



enhance cloud-cloud collision, thus increase both star for- 
mation and accretion in a self-gravitation disk. According 
to the numerical results of Wang & Biermann (1998) for 
the QSO evolution, we adopt the assumption (3x^2 ^ 1 
(i.e. the star formation approximately scales with the ac- 
cretion). In this case, (32 ^ 20, corresponding to a star 
formation rate WMq/yt in the central disk, which is 
at a reasonable level for the observed star formation rate 
in Seyferts. The initial seed black hole in our model is tiny, 
~ 5 Mq. We know from the calculation that the black hole 
grows by a rate oc M^j^ at the beginning till i?acc ^ Rm, 
afterwards, the accretion rate would decrease by 1/Afbh 
and only a mediate black hole of ~ 10^ Mq can be formed 
if there is no external trigger. In this case, a tidal pertur- 
bation around 5 lO^yr could help to grow a massive black 
hole, which corresponds to z 0.8 in an Einstein-de Sitter 
Universe with fl = I. Hq = 50kms^^ Mpc~^ (where we 
think the tidal interactions happen frequently, also could 
be a peak of Seyfert activities) . 

Although the parameters, such as initial black hole 
mass, the inner boundary of the accretion disk etc. could 



have certain influence on the mass ratio, we notice from 
the numerical results and the illustrative estimation that 
velocity dispersion of the accreting gas is another key pa- 
rameter besides the value of /?i/?2, which could cause a 
large dispersion of the black hole to bulge mass ratio in 
AGNs. Fig. shows a correlation between the nuclear 
black hole to bulge mass ratio and the velocity disper- 
sion of accreting gas to the center from our simulation, 
which gives a rough proportionality of Afbh/Mbuige oc vl '^. 
In our model, we are aiming to explain a broad distribu- 
tion of black hole to bulge mass ratio of a similar bulge 
system in QSOs and Seyferts. If we adopt the masses of 
the spheroidal systems scale with their luminosities as 
M oc L^/'* and the scale length r oc L^/^, combining 
with the virial theorem gives the bulge velocity disper- 



sion oc 4f oc oc oc M^/^ (Faber & Jackson 



1976, Faber et al. 1987, Peterson 1997). So, we obtain 
.^^^buigo OC a^'^. The assumption that the velocity disper- 
sion of accreting gas is systematically proportional to the 
bulge velocity dispersion gives a correlation of Mbh oc w^'^, 
which is shallower than the slope of Mbh oc a^'^ given by 
Ferrarese & Merritt (2000), but steeper than the best-fit 
correlation (Mbh oc cr^ "^^) by Gebhardt et al. (2000). We 
see from Fig. |^ that the variation of the velocity dispersion 
of accreting gas would cause a broad distribution of black 
hole to bulge mass correlation, where the accreting gas 
with higher velocity dispersion close to the bulge veloc- 
ity dispersion could lead to an accretion near Eddington 
limit, and form more massive black holes in QSOs than 
in Seyferts of similar bulges. Therefore, the black hole 
to bulge mass correlation in massive systems is close to 
an upper limit and much tighter than that in medium- 
bulges. The physical interpretation of such discrepancy 
could be due to the different formation or evolution envi- 
ronment of the two kinds of systems (Seyferts and QSOs). 
The violent collision between two galaxies could trigger 
intense starburst in the center, heat or shock the interstel- 
lar medium efficiently than tidal interaction, thus drive a 
Bondi flow fuelling the central black hole with probably a 
higher sound speed and result in a higher accretion rate. 
In this case, it may enhance a QSO evolution with the 
accretion rate close to Eddington limit, form a massive 
black hole, and lead to a black hole to bulge mass ratio 
higher than in case of Seyferts. 

3. Summary 

Recent reverberation mapping of a sample of Seyfert 1 
galaxies by Wandel (1999) suggests that the black hole to 
bulge mass correlation in Seyferts has a significant disper- 
sion with a mean value of ~ 10~^'^, which is about one 
magnitude lower than in case of QSOs. Considering a sim- 
ple unified formation scheme for AGNs, we demonstrated 
a possible black hole evolution scheme where we assume 
black hole accretes gas coming within its influence by a 
uniform Bondi flow. This scenario could interpret not only 
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the statistical mass correlation in AGNs and normal galax- 
ies, but a large dispersion of the black hole to bulge mass 
ratio in QSOs and Seyferts. We found the black hole to 
bulge mass correlation in such evolution scheme strongly 
depends on the velocity dispersion of the accreting gas, 
which might be an important environmental parameter 
for the observed correlation besides the relation between 
the star formation and accretion in the central accretion 
disk during galaxy interactions (Wang & Biermann 1998, 
2000). Thus, the resuhs of Wang & Biermann (1998) rep- 
resent a limiting case for the black hole evolution in AGN, 
where the accretion is close to Eddington rate. The black 
hole to bulge mass ratio versus velocity dispersion of ac- 
creting gas to the central region is plotted in Fig. |l|, which 
gives a rough proportionality of Mbh/Afbuigo v^-* for a 
fixed bulge system in our model. If we consider a relation 
of bulge mass with the galaxy velocity dispersion in the 
virial equilibrium (Afbuigo oc a^'^), we could estimate a 
correlation of the nuclear black hole mass with the veloc- 
ity dispersion of the accreting gas as Mhh oc w^'^, within 
the slope suggested by recent work of Ferrarese & Merritt 
(2000) and Gebhardt et al. (2000). The square and the 
star in Fig. ^ correspond to the mean value of the black 
hole to bulge mass ratio in two typical systems ( Seyferts 
and QSOs), which shows the velocity dispersion of the 
accreting gas to form a QSO phase is much higher than 
in case of a Seyfert, almost close to the virial velocity of 
the system. We present a possible interpretation for such 
scenario, where we think the intense starburst during vio- 
lent mergers could heat or shock the interstellar medium 
(ISM) to a higher sound speed more efficiently than by 
tidal interactions, result in a higher accretion rate close 
to Eddington limit and a higher black hole to bulge mass 
ratio finally in QSOs and early type galaxies than in case 
of Seyferts. 
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